Variability of the H/5 line profiles as an indicator of orbiting bright 
spots in accretion disks of quasars: a case study of 3C 390.3 
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ABSTRACT 

Here we show that in the case when double peaked emission hnes originate from outer parts 
of accretion disk, their variabihty could be caused by perturbations in the disk emissivity. In 
order to test this hypothesis, we introduced a model of disk perturbing region in the form of a 
single bright spot (or flare) by a modification of the power law disk emissivity in appropriate 
way. The disk emission was then analyzed using numerical simulations based on ray-tracing 
method in Kerr metric and the corresponding simulated line profiles were obtained. We applied 
this model to the observed H/3 line profiles of 3C 390.3 (observed in the period 1995-1999), 
and estimated the parameters of both, accretion disk and perturbing region. Our results show 
that two large amplitude outbursts of the H/3 line observed in 3C 390.3 could be explained by 
successive occurrences of two bright spots on approaching side of the disk. These bright spots 
are either moving, originating in the inner regions of the disk and spiralling outwards by crossing 
small distances during the period of several years, or stationary. In both cases, their widths 
increase with time, indicating that they most likely decay. 
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The huge amount of Active Galactic Nucleus 
(AGN) energy is released through accretion onto 
super-massive black hole (BH), supposed to exist 
in the center of AGN. The emission of the accre- 
tion disk is not only in the continuum, but also in 
the emission lines (e.g. in Fe Ka line) and in low 
ionization lines, as e.g. in broad Balmer emission 
lines which are seen as double peaked (DP). DP 
Balmer lines are found in 20% of radio loud AGN 



at z < 0.4 (|Eracleous fc Halpern lHool |2003| ) and 
4% of the Sl oan digital Sky Survey (SDSS) quasars 



20031) . 



at z < 0.33 (jStrateva et al. 

Broad, double-peaked emission lines of AGN 
provide dynamical evidence for presence of an ac- 
cretion disk feeding a supermassive black hole in 
the center of AGN. But in some cases, the variabil- 
ity of these lines shows certain irregularities which 
could not be explained just by standard model of 
accretion disk. 

The DP line profiles are often used to extract 



the disk parameters (see e.g. Chen et al. _ 1989t 



Chen fc Halpern llQQd 'Eraclcous fc Halpern Ill994 . 
20031 : IStrateva et al. 12 008: Eraclco us et al. II2009I ). 



In a series of papers Dumont fc CoUin-Souffrin (see 
CoUin-Souffrin ' 'igSV"; 'Collin-Souffrin fc Du moirt" 
1990; Dumont & Collin-Souffrin 1990a, be, and 
references therein) investigated the radial struc- 
ture and emission of the outer regions of the opti- 
cally thin accretion disks in AGN and calculated 
detailed grid of photoionisation models in order 
to predict the relative strengths of low-ionization 
lines emitted from the disk. They obtained inte- 
grated line intensities and line profiles emitted at 
each radius of the disk, for its different physical 
parameters. They also studied the influence of 
the external illumination on the structure of the 
disk, considering the point source model, where 
a compact source of non-thermal radiation lo- 
cated at a given height illuminates the disk and 
the diffusion model, where the radiation of a cen- 
tral source is scattered back towards the disk by 
a hot diffusing medium. One of the first meth- 
ods for calculating the profiles of optical emission 
lines fro m a relativist i c accr etion disk was pro- 
posed bv IChen et al. I ( 19891 ). The limitation to 
this method is that the accretion disk structure 
required to explain the variability of the line pro- 
files cannot be axi-symmetric, i.e. very often the 
red peak is higher than blue one and that can- 
not be explained by this model. It is not possible 
in a circular disk, in which the blue peak is al- 
ways Doppler bo osted to be stronger tha n the red 
peak. Therefore, lEracleous et al. I (|l995l ) adapted 
the circular accretion disk model to elliptical disks 
in order to fit the profiles of double-peaked emit- 
ters with a red peak stronger than the blue one. 
This model introduced eccentricity and phase an- 
gle parameters to the circular model described 
above, and the pericenter distanc e of the elliptical 



orbits (see lEracleous et al. 1119951 ) 



Spectroscopic monitor ing of double-peaked 

emitters (see e.g.lShapovalova et al. 2001; Goza ri et al. 
200?! : IShapovalova et al. 1120091) has revealed that 
a ubiquitous property of the double-peaked broad 
emission lines is variability of their profile shapes 
on the timescales of months to years. DP line pro- 
files are observed to vary on timescales of months 
to years, i.e. on timescales of the order of the 



dynamical time or longer (e.g. IVeilleux fc Zheng 
19911 : I Zheng et all Il99ll : iMarziani et al. I Il993t 



Roniano et al. ll998l:ISergeev et al. 2000l:IShapova,lova et al 



2001 : IStorchi-Bergmann et al. ll2003l : lGezari et al 



2003). 

This slow, systematic variability of the line pro- 
file is on the timescale of dynamical changes in an 
accretion disk, and has been shown to be unre- 
lated to the shorter timescale variability seen in 
the overall flux in the line, due to reverberation 
of the variable ionizing continuum. Patterns in 
the variability of the broad Balmer lines are of- 
ten a gradual change and reversal of the relative 
strength s of the blue-shifted an d red-shifted peaks 
(see e.g. lNewman et al. Ill997l ). 

Periodic variability of the red and blue peak 
strengths has also been attributed to a precess- 
ing elliptical disk, a precess ing single-armed spira l 
(as e.g. 3C 332, 3C 390.3: iGilbert et a l. ' ('1999'); 
NGC 1097: IStorchi-Bergmann et al.T^OOa) ). and 
a precessing wa rp in the disk. For instance, 
Wu et al. I ( 20081) computed the proflles of Balmer 
emission lines from a relativistic, warped accre- 
tion disk in order to explore the certain asymme- 
tries in the double-peaked emission line profiles 
which cannot be explained by a circular Keplerian 
disk. Elliptical disks and spiral w aves have been 



detec t ed in cataclysmic variables (jSteeghs et al. 



199?! : iBaptista fc Catalan I |2000|), and a radia 



tion induced warp has been detec ted in the large- 



scale disk of the AGN NGC 4258 (|Malonev et al. 
19961) . 



Spiral waves are a physically desirable model 
since they can be produced by instability in the 
vicinity of a black hole. They can play an impor- 
tant role in accretion disks, because they provide 
a mechanism for transporting angular momentum 
outward in the disk, allowing the gas to flow in- 
wards, towards the central black hole. Long-term 
profile variability is thus a useful tool for extract- 
ing information about the structure and dynam- 
ics of the accretion disk most likely producing the 
double-peaked emission lines. 

In this paper, we present an investigation of 
the disk line variations due to instability in accre- 
tion disk. First we developed a model, assuming 
that instability in the accretion disk affects disk 
emissivity. This model and some simulations of 
expected line profile variability are presented in 
§2. In §3 we compare the model with observa- 
tions taken from long-term monitoring of 3C 390.3 



tions taken trom long-terrn monitormg oi d(J dyu.i 
( Shapovalova et al. 1120011 ) in order to obtain pa- 



rameters of perturbations. In §4 we discuss our re- 
sults in the light of possible physical mechanisms 
which could cause such perturbations, and finally, 
in §5 we outline our conclusions. 

2. The model of perturbation in the accre- 
tion disk 

Here we introduce the model and some approx- 
imations used in the simulations of accretion disk 
perturbation. 

2.1. Long term variation of DP line pro- 
files: some assumptions and problems 

As we mentioned above, the DP line profile vari- 
ability does not appear to correlate with changes 
in the line and/or continuum flux, and conse- 
quently one can assume that the changes in line 
profile are likely caused by changes in the accre- 
tion disk structure. There are several examples 
of the long-term variability (on timescales of sev- 
eral years) of the DP line profile of some objects 
which has been successfully modeled by the pre- 
cession of a non-axisymmetric accretion disk, such 
as an elliptica, l disk o r a disk with a spiral arm 



(iGczari et aTl l2007t 



Storchi-Ber gmann et al. 
2003; S hapovalova et al. I I2OOII: ,Gilbcrt et al. 



19991 and references therein). These models, how- 
ever, fail to explain the long-term variability of 
some objects and the short-term variability (on 
timescal es from seve r al months to a ye ar) of all 
objects JLcwis "2005'; 'Lew is et al."1l2010l) . For in- 
stance, [Lewis et al. (2010() found that the two 
simple models, an elliptical accretion disk and 
a circular disk with a spiral arm, are unable to 
reproduce all aspects of the observed variability, 
although both account for some of the observed 
behaviors. Therefore, these authors suggest that 
many of the observed variability patterns could 
be reproduced assuming a disk with one or more 
fragmented spiral arms. 

Other attempts to explain the DP line profile 
variability through perturbations of the disk struc- 
ture introduced bright spots ove r an axisymmet 



ric acc retion disk. As an example iNewman et al. 



(|l997[ ) successfully modeled the variation of the 
Ha peak intensity ratio of Arp 102 B with a single 



spot rotating within the disk, but iGezari et al 



In the case of Fe Ka variability [Turner et al. 
pood ) used the spot model to explain the vari- 
ability of the iron li ne profile of M rk 766 in the 
X-ray band. Also, iDovciak et al. I (|2008[ ) stud- 
ied variations of the iron line due to an orbit- 
ing spot which arise by reflection on the surface 
of an accretion disk, following its illumination by 
an X-ray fiare in form of an off-axis point-like 
source just above the accretion disk. Besides the 
spots in accretion disk, the Fe Ka line of some 
AGN could be also significantly affected by highly- 
ionized fast accretion disk outfiows. For instance, 
Sim et al. I ( 20101 ) found that the major features 
in the observed 2-10 keV spectrum of the bright 
quasar PG1211-I-143 can be well reproduced by 
their Monte Carlo radiative transfer simulations 
which include a variety of disk wind (outfiow) 
models. 

To explain th e short-timescale vari a bility of the 
DP line profiles iFlohic fc Eracleousi ( 2008[ ) con- 



structed stochastically perturbed accretion disk 
models and calculated Ha line profile series as the 
bright spots rotate, shear and decay. They ruled 
out spot production by star/disk collisions and fa- 
vor a scenario where the radius of marginal self- 
gravity is within the line emitting region, creating 
a sharp increase in the radial spot distribution in 
the outer parts. 

2.2. The model of bright spot like per- 
turbing region 

We model the emission from accretion disk us- 
ing numerical simulations ba sed on ray-tracing 
method in Kerr metric (see e.g. lJovanovic fc Popovic 
2009ai and references therein). Although this 
method was developed for studying the X-ray 
radiation which originates from the inner parts 
of t he disk close to the central black hole (see 
e.g. Ijovanovic fc Popovic I l2008al ). it can be also 
successfully applied for modelling the UV/optical 
emission which originates from the outer regions 
of the dislfl. 



(|2007[ ) were not able to apply the same model to 
the same object at a different time period. 



^ As shown in Ijovanovic fc Popovic I ll2008al ) , the effects of 
strong gravitational field and angular momentum of rotat- 
ing black hole are significant only in the innermost parts 
of accretion disk, in vicinity of the central supermassive 
black hole, up to the several dozens of gravitational radii. 
In the outer parts of the disk, such as those where H/3 line 
originates, these effects are negligible and Kerr metric with 
zero angular momentum, i.e. Schwarzschild metric, is a 
very good approximation. 



Surface emissivity of the disk is usually as- 
sumed to vary with radius as a power law (e.g. 



Popovic et al. 1120031 ): e(r) = Eq ■ r*, where Eq is 



an emissivity constant and q - emissivity index. 
Total observed flux is then given by: 



FobsiEobs)= I e{r) 

image 



/e 



Eobs - gEo 



(1) 



where g is the energy shift due to the relativistic 
effects: g — — ^^, i?o is the rest energy of the line, 

a is the local turbulent broadening and d^ is the 
solid angle subtended by the disk in the observer's 
sky. 

In this paper we adopt the following modifi- 
cation of the power-law disk emissivity, in order 

to intmrliirp a hrigtif gpnt.-1ilrp pprtiirhirig T-pgij 



in the disk (jJovanovic fc Popovic ll2008Ul2009 



Stalevski et al. 1 120081 ): 



m 



£\{x,y) =£{r{x,y)) 
■ \ + £p ■ e 



y-Vv 



(2) 

where ei{x,y) is the modified disk emissivity at 
the given position {x, y) expressed in gravitational 
radii Rg, e{r{x,y)) is the ordinary power-law disk 
emissivity at the same position, Sp is emissivity 
of the perturbing region (i.e. amplitude of the 
bright spot), (xp,yp) is the position of perturb- 
ing region with respect to the disk center (in Rg) 
and {wx,Wy) are its widths (also in Rg). A 3D 
plot of above expression for modified emissivity 
law is given in Fig. [T] This simple model is 
suitable for our purpose because it allows us to 




Fig. 1. — A 3D plot of modified disk emissivity 
given by Eq. (0) for 200 Rg < r{x, y) < 1200 Rg, 
q = —2.5 and for the following parameters of 
perturbing region: Sp = 5, Xp = 700 Rg, yp = 



change amplitude, width and location of bright 
spot in respect to the disk center. In that way we 
are able to simulate displacement of bright spot 
along the disk, its widening and amplitude de- 
crease with time (decay). Moreover, the above 
bright spot model can be successfully applied for 
studying the variations of accretion disk emission 
in different s pectr al band s, from X- rays to optical 
band (see e.g.ljoya novic fc Popovic ..2008b . ,2009bt 



Stalevski etall 120081 ). 
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100 Re 



3. Results: Model vs. observations 

3.1. Perturbation in the accretion disk: 
modeled profiles 

In order to test how this model of bright spot 
affects the H/3 line profile, we performed several 
numerical simulations of perturbed emission of an 
accretion disk in Kerr metric for different posi- 
tions of bright spot along x and y-axes in both, 
positive and negative directions. For these sim- 
ulations we adopted the following parameters of 
the disk: inclination i = 30°, inner and outer radii 
Rin — 200 and Rout = 1200 Rg, power law emis- 
sivity with index g = — 1, local turbulent broaden- 
ing a = 2000 km s~^ and normalized angular mo- 
mentum of black hole a = 0.5. The corresponding 
results are presented in Fig. [21 As it can be seen 
from this figure, when the bright spot moves along 
the positive direction of x-axis (receding side of the 
disk) it affects only " red" wing of the line (Fig. [21 
top right), but when it moves along the negative 
direction of x-axis (approaching side of the disk) 
it affects only "blue" wing of the line (Fig. [21 
top left). In both cases, the other wing and the 
line core stay nearly constant, and therefore al- 
most unaffected by bright spot. The situation is 
quite opposite when the bright spot moves in both 
directions along the y-axis, because then it affects 
only the line core, while the both of its wings stay 
almost intact (see the bottom panels of Fig. [2]) . 

We also performed the corresponding simula- 
tions for different positions of a bright spot which 
moves from the inner radius of the accretion disk 
towards its outer parts along the y = x direction, 
and found similar behavior of the simulated line 
profiles (see Fig. [3]). As one can see from Fig. 
[3l for certain positions of perturbing region along 
y — X direction we obtained the line profile with 
almost symmetrical wings, while in other cases ei- 
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Fig. 2. — Shapes of perturbed emissivity of an accretion disk in Kerr metric and the corresponding perturbed 
(dashed hue) and unperturbed (sohd hne) H;5 hne profiles for the following parameters of perturbing region: 



-p 



Uy 



200 Rg. The positions of perturbing region are: 



Xp = 700 Rg and yp 
(bottom right). 



= -700 Rg and Vp = (top left), 
(top right), Xp — and yp — —700 Rg (bottom left) and Xp — and yp — 700 Rg 
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Fig. 3. — Variations of the perturbed H/3 line profile for different positions of bright spot along the y = x 
direction. Disk parameters are the same as used in examples in Fig. [2] Left panel corresponds to the 
positions of bright spot on approaching side of the disk, while the right panel corresponds to the receding 
side of the disk. In both cases, the positions of bright spot were varied from the inner radius (bottom profiles) 
of the disk towards its outer radius (top profiles) . 



ther the "blue" peak is brighter than the "red" 
one, or the "red" peak is stronger than the "blue" 
one. 

The next step in our analysis was to use our nu- 
merical simulations for fitting the observed spectra 
of 3C 390.3 in order to study the variability of its 
H/3 spectral line due to emissivity perturbations 
in its accretion disk. 

3.2. Observations of 3C 390.3 

To test the model, we used 22 spectra of 3C 
390.3 observed from Novemb er 1995 until June 
1999 (see Fig. 8 in [Shapovalov a et al. | [2001'). 

Spectra of 3C 390.3 were taken with the 6 
m and 1 m telescopes of the SAO RAS (Rus- 
sia, 1995-2001) and at INAOE's 2.1 m telescope 
of the "Guillermo Haro Observatory" (GHO) at 
Cananea, Sonora, Mexico (1998-1999) in moni- 
toring regime in 1995-1999. They were obtained 
with long slit spectrographs, equipped with CCD 
detector arrays. The typical wavelength interval 
covered was that from 4000 A to 7500 A, the spec- 
tral resolution varied between 5 and 15 A and the 
S/N ratio was > 50 in the continuum near Ha and 
H/3. Spectrophotometric standard stars were ob- 
served every night. The spectrophotometric data 
reduction was carried out either with software de- 
veloped at SAO RAS or with the IRAF package for 
the spectra obtained in Mexico. The image reduc- 
tion process included bias, flat-field corrections, 
cosmic ray removal, 2D wavelength linearization, 
sky spectrum subtraction, addition of the spec- 
tra for every night, and relative fiux calibration 
based on standard star observations. Spectra were 
scaled by the [O III] AA4959+5007 integrated hue 
flux under the assumption that the latter did not 
change during the time interval covered by our 
observations ( 1995-2001) . A value of 1.7x10"" 



ergs s ^cm ^ (jVeilleux fc Zheng II 199 II ) for the in- 



tegrated [O III] line flux was adopted. In order 
to calculate a normalization coefficient, the con- 
tinuum was determined in two 30 A wide clean - 
line free- windows centered at 4800 A and 5420 A. 
After continuum subtraction, blend separation of 
the H/3 and [O III] components was carried out by 
means of a Gaussian fitting procedure, applied to 
the following: H/3 — broad blue, broad red and 
central narrow; [O III] AA4959,5007 — broad and 
narrow components. The forbidden lines are rep- 
resented by two Gaussian curves with an intensity 



ratio I(5007)/I(4959)=2.96. 

Comparisons between mean and rms spectra of 
the Ha and H/3 broad line profiles of 3C 390.3 dur- 
ing 1995-2007 (including the H/ 3 spectra from this 
paper ) are given in Fig. 12 of 'Shapo valova et al 
( 20101 ). It can be seen that these profiles are sim- 



ilar, and moreover , the corresponding Ha pr ofiles 



of 3C 390.3 from Shapovalova et al.l ( 20101 ) and 



Gezari et al. I ( 20071 ) are also similar. These com- 



parisons indicate that the [O III] AA4959,5007 nar- 
row lines subtraction in H/3 region was performed 
correctly. 

The flux of H/3 and the broad component profile 
was obtained from scaled spectra after continuum 
subtraction and removal of the [O III] doublet and 
narrow H;9 component. Then the observed con- 
tinuum fluxes were c orrected for the apertu re ef- 
fects using scheme by |PetersoiLet_alJ (1995|). The 
mean error (uncertainty) in our flux determina- 
tions for the continuum and H/ 3 flux is <3%. More 



details can be found in paper (Shapovalova et al. 
(I2OOII) . 



3.3. 



Perturbation in the accretion disk of 
3C 390.3 



Radial velocities of the blue and red peaks 
of the H/3 and Ha broad lines of 3C 390.3 vary 
with time (see e.g.iGaskell Ill996t Eracleous et al 



1997 : Shapovalova et al.l2010r ). IShapovalova et al 
(|200 1l obtained the H/3 difference profiles by sub- 
tracting the average spectrum corresponding to 
the minimum activity state (September 9, 1997) 
from the individual spectra (see their Fig. 11). 
These authors found (see their Table 7) that the 
radial velocity of the blue peak increased from 
-3200 km s-i in 1995-1996 to -5200 km s"^ in 
1999. At the same time the radial velocity of the 
red peak increased from -1-4900 km s"^ in 1995- 
1996 to -1-7000 km s'^ in 1999. Here we analyze 
the possibility that the velocities, corresponding 
to the peak shifts in the H/3 integral and differ- 
ence profiles, vary with time due to perturbations 
in disk emissivity. 

In order to fit the spectral H/3 line shapes of 3C 
390.3, we first estimated the disk parameters from 
several profiles. We found following parameters of 
the disk: inclination i = 20°, inner and outer radii 
Rin = 100 and Rout = 1300 Rg, broken power law 
emissivity with index q = —\ for Rin < r < R^r 



and q 



-3 for Rhr < r < Rout, radius at which 



= 500 Eg, emis- 
1, local turbulent 



slope of eniissivity changes Rbr 
sivity of perturbing region Sp — 
broadening a = 2000 km s~^ and normalized an- 
gular momentum of black hole a — 0.5. These val- 
ues of parameters are in accordance with the cor- 
responding parameters for 3C 390.3 obtained by 
Flohic fc Eracleous I ()2008l ). who found the follow- 
ing values: i = 27°, i?„ = 450, Rout = 1400 Rg 
and a = 1300 km s~^. In all our simulations 
we held emissivity of perturbing region fixed at 
Ep — 1, but nevertheless, the brightness of per- 
turbing region varies with time due to its change- 
able area (i.e. its widths are taken as free pa- 
rameters which vary with time), as well as due to 
power law emissivity of the disk which decreases 
with radius. We studied stationary and moving 
perturbing regions, and in the latter case their co- 
ordinates are also assumed as free parameters. 

The fitting of the observed H/3 line shapes of 3C 
390.3 is performed separately in the case of moving 
and stationary perturbing regions in the following 
three steps: (i) the simulated H/3 line profiles are 
calculated for some initial set of values of free pa- 
rameters; (ii) the corresponding root mean square 
(RMS) of residuals between the observed and sim- 
ulated line profiles is calculated; (iii) the values of 
free parameters are varied and the procedure (i)- 
(iii) is repeated until the RMS deviations (RMSD) 
between the observed and fitted H/3 profiles be- 
come as small as possible. 

Fig. m shows comparisons between all 22 ob- 
served spectra (black solid line) and the obtained 
best fits assuming the moving (red solid line) and 
stationary (blue dashed line) perturbing regions, 
while the corresponding positions of the perturb- 
ing regions are presented in Fig. [5j As it can 
be seen from Fig. \5\ the obtained positions of 
the moving perturbing regions are distributed in 
the form of two spiral arms, indicating that these 
perturbations originate in the inner regions of the 
disk and spiral away towards its outer parts, mov- 
ing faster in the azimuthal direction as they get 
further away from the center of the disk. Table [T] 
contains the obtained values of the fitted param- 
eters, as well as the corresponding RMSD in the 
case of both, moving and stationary perturbing re- 
gions. We first studied the moving perturbations, 
but after examining the obtained results we found 
that, during a period of a few years, such perturba- 



tions move only by small distances which are com- 
parable to their widths (see Tabled]). Therefore, 
in order to test whether the obtained displace- 
ments are reliable, we repeated the fitting, but this 
time assuming the stationary perturbing regions 
with variable widths. In this case we obtained 
the best fits for two perturbations positioned at 



= -100 R 



•9' 



y 



220 Rg and at x 



-220 R 



ff' 



y ~ 125 Rg, respectively (denoted by crosses in 
Fig. [5]). As one can see from Fig. |31 and Table 
[1] regardless the significant variations of the H/3 
line profile during the analyzed period, both mod- 
els of perturbing region resulted in similar fits for 
the most of these spectra, except for the spectra 
observed during 1999, where the moving perturb- 
ing region achieved better fits. The latter result 
should be taken with caution because the obser- 
vations from 1999 were performed after a large 
gap of « 1 year. Therefore, the displacements of 
the perturbing regions cannot be considered as in- 
disputably confirmed, but on the contrary, their 
widths almost certainly vary with time. The last 
conclusion is valid for both, moving and station- 
ary perturbing regions, since neither of them can 
provide satisfactory fits with fixed widths. 

The obtained best fit positions of both, moving 
and stationary perturbations are located on the 
approaching side of the accretion disk (see Fig. 
[5|) and these perturbations can be most likely at- 
tributed to successive occurrences of two different 
bright spots. This assumption is in good agree- 
ment with observations, since two large ampli- 
tude outbursts of H/? line are observed duri ng the 
analyzed period (jShapovalova et al. 1 120011 ) , and 
therefore each bright spot can be assigned to one 
of them: the bright spot which positions are de- 
noted by squares corresponds to October 1994 - 
July 1997 outburst while the other one, which po- 
sitions are denoted by triangles, corresponds to 
July 1997 - June 1999 outburst. Using the time 
differences between two successive observed spec- 
tra we were able to estimate the speeds of both 
moving bright spots (see Table [ij. For an average 
velocity of the first bright spot we obtained the 
value of 7298 km s^^ and for the second one 6575 
km s^^. As it can be seen from Table [1] widths of 
bright spots are increasing with time, indicating 
that they decay until they completely disappear. 
It should be noted that, inevitably, there is a cer- 
tain degree of degeneracy in the parameter space. 
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Fig. 4. — Comparisons between the observed H/3 line profiles of quasar 3C 390.3 (black solid line) and 
the corresponding simulated profiles due to two successive bright spots. The red solid line represents the 
simulated profiles due to the moving bright spots which estimated positions are presented in Fig. [5] and 
widths in Table [T] The blue dashed line corresponds to the simulated profiles due to the stationary bright 
spots, positioned at a; = —100 Rg, y = 220 Rg during the first outburst and at a; = —220 Rg, y = 125 Rg 
during the second outburst. 
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Fig. 4. — (continued) 
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Fig. 5. — Positions of moving perturbing region along the accretion disk corresponding to two observed 
amplitude outbursts: October 1994 - July 1997 outburst (squares) and July 1997 - June 1999 outburst 
(triangles). In both cases, the moving perturbation originates in the inner regions of the disk and spirals 
away towards its outer parts. For average speed of the first perturbation we obtained value of 7298 km s~^ 
and for the second one 6575 km s~^. Positions of the stationary perturbing region (denoted by crosses in 
the right panel) are: x — —100 i?g, y — 220 Rg during the first outburst and x — —220 Rg, y — 125 Rg 
during the second outburst. 



since similar results could be obtained with some- 
what different combinations of perturbing region 
positions and widths. 

4. Discussion 

Several physical mechanisms could be responsi- 
ble for perturbations in accretion disk emissivity, 
i.e. for bright spot formations. The most plau- 
sible candidates for such mechanisms are: disk 
self-gr avity, baroclinic vorticity, disk-star colli- 



sions (|Flohic &: Eracleous I l20Q8l and references 



therei n), tidal disruptions of stars by central black 
hole (jStrubbe fc Quataert I l2009l a nd references 



there in) and fragmented spiral arms ([Lewis et al. 
I2OIOI and references therein). 

The disk self-gravity, driven by Jeans instabil- 
ity, could cause production of clumps in the disk 
which have typical sizes in the range from 10 to 
1000 gravitational radii for a 10^ M© central black 
hole. Such clumps do not shear with differential 
rotation and they have high brightness that varies 
very little over time. Since the obtained results 
indicate that bright spots decay by time and spi- 
ral along the disk, it is not likely that these bright 
spots could be identified as clumps created by self- 
gravity, although their sizes are comparable. 

Baroclinic vorticity appears in the accretion 



disk due to its differential rotation in combination 
with the radial temperature gradient, causing the 
material in the disk to spiral around the center of 
vortex. Such a vortex would have higher density, 
and hence higher brightness, causing the forma- 
tion of a bright spot. The typical sizes of such 
spots, as well as their shearing with differential 
rotation of the disk, are still unknown since differ- 
ent numerical simulations ga ve contradictory re- 



sults (for more details see e.g. iFlohic fc Eracleous 



20081 and references therein). Therefore, in the 



case of 3C 390.3, this mechanism still cannot be 
neither ruled out, nor accepted. 

Disk-star collisions are assumed to be very fre- 
quent events which happen on daily timescales and 
which could increase disk surface temperature in 
the region of collision, and thus, create a bright 
spot. Such bright spots shear with differential ro- 
tation of the disk and decay as the material cools 
down. However, the typical size of a such bright 
spot immediately after collision is close to the size 
of the star, which is very small when expressed in 
gravitational radii. Therefore, neither this mech- 
anism could be accepted as a potential cause of 
two bright spots, detected in the case of quasar 
3C 390.3. 

Th e tidal disruption of stars by central black 
hole (jStrubbe fc Quataert I l2009l and references 
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Table 1 
Parameters of perturbing region obtained by fitting the observed spectra. 



Date 


JD (2400000+) 


X[Rg] 


y[R.] 


W[Rg] 


d[R,] 


w[km s "*"] 


RMSD 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


The First Outburst 


















1995 Nov 17 


50039.156 


-30 


175 


100 


100 






0.09730 


0.09460 


1996 Feb 14 


50127.602 


-40 


200 


105 


130 


26.93 


5203.24 


0.09680 


0.09900 


1996 Mar 20 


50162.580 


-50 


208 


106.5 


140 


12.81 


6257.63 


0.08476 


0.09138 


1996 Jul 12 


50276.567 


-75 


230 


110 


180 


33.3 


4993.36 


0.08518 


0.08178 


1996 Jul 17 


50281.434 


-76.3 


231.3 


110 


180 


1.84 


6456.24 


0.08820 


0.08249 


1996 Aug 10 


50305.489 


-90 


235 


112 


190 


14.19 


10082.9 


0.09435 


0.09201 


1996 Sep 11 


50338.309 


-115 


238 


120.5 


200 


25.18 


13112.6 


0.12337 


0.12426 


1997 Mar 04 


50511.622 


-165 


245 


130 


260 


50.49 


4978.94 


0.11448 


0.12082 


The Second Outburst 


















1997 Aug 30 


50691.463 


-160 


100 


100 


100 






0.13717 


0.13593 


1997 Sep 10 


50701.576 


-165 


105 


105 


110 


7.07 


11950.55 


0.12034 


0.11905 


1997 Dec 30 


50813.195 


-205 


120 


145 


175 


42.72 


6541.48 


0.09499 


0.09398 


1998 Jan 22 


50835.631 


-215 


122.5 


155 


195 


10.31 


7852.39 


0.10690 


0.10949 


1998 Feb 23 


50867.560 


-225 


125 


165 


215 


10.31 


5517.75 


0.09825 


0.09618 


1998 May 06 


50940.354 


-245 


129 


185 


265 


20.40 


4788.88 


0.08867 


0.08965 


1998 Jun 25 


50990.302 


-265 


132 


205 


285 


20.22 


6920.32 


0.07229 


0.07387 


1998 Jul 16 


51010.719 


-275 


134 


210 


295 


10.20 


8537.04 


0.06809 


0.07141 


1998 Jul 25 


51019.723 


-278 


134.5 


212 


300 


3.04 


5773.22 


0.06548 


0.06615 


1998 Aug 30 


51055.551 


-285 


135.5 


215 


320 


7.07 


3373.22 


0.07590 


0.07638 


1998 Sep 26 


51082.429 


-293 


136 


218 


340 


8.01 


5097.09 


0.07948 


0.07788 


1999 Aug 19 


51410.309 


-500 


140 


250 


400 


207.04 


10792.43 


0.10645 


0.11163 


1999 Sep 04 


51426.208 


-505 


140 


252 


400 


5 


5375.06 


0.09147 


0.09232 


1999 Oct 03 


51455.172 


-510 


140 


254 


400 


5 


2950.49 


0.09977 


0.10270 



Note. — Col. (1): date of observation; Col. (2): epoch of observation (in JD); Cols. (3-4): x and y 
coordinates of the moving bright spots; Col. (5): widths (w = Wx = Wy) of the moving bright spots; Col. (6): 
widths {w = Wx ~ Wy) of the stationary bright spots, positioned at a; = —100 Rg, y — 220 Rg during the first 
outburst and at a; = —220 Rg, y — 125 Rg during the second outburst; Col. (7): linear distance crossed by 
the moving bright spots between two successive observations; Col. (8): average speed of the moving bright 
spots between two successive observations; Col. (9); root mean square deviation between the observed and 
fitted H/3 line profiles in the case of the moving bright spots; Col. (10): the same as Col. (9), but in the case 
of the stationary bright spots. 
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therein) happens when a star passes the tidal ra- 
dius of the black hole, i.e. when the black hole's 
tidal gravity exceeds the star's self-gravity. Gas 
of a disrupted star falls back to the black hole 
at a super-Eddington rate, releasing a flare of en- 
ergy which then blows away a significant frac- 
tion of the falling gas as an outflow. Such super- 
Eddington flares and outflows could induce insta- 
bilities in accretion disk in form of bright spots. 
However, this mechanism is an unlikely candi- 
date for a potential cause of the bright spots in 
the case of 3C 390.3 due to the following weak- 
nesses: (i) super-Eddington outflows are short liv- 
ing (^ 10 days); (ii) frequency of star disrup- 
tions in a typical elliptical galaxy is very low, be- 
tween 10^^ and 10^^ per year, and in the case 
of the black hole of 3C 390.3 w hich mass is ^ 
5 X IO^Mq (JLewis fc Eracleous ||2006). it i s near 



the low end of this range (see e .g. 



Wang fc Merritt 



2004 iMagorrian fc Tremaine Ill999t ): (iii) anv stel- 
lar debris could be hardly released in the case of 
3C 390.3 since the main sequence stars are dis- 
rupted within the innermost stable orbit around a 
non-rotating black hole as massive as that in 3C 
390.3; (iv) it is unclear how the bright spots pro- 
duced by this mechanism should evolve with time. 
Spiral arms in AGN disks could be formed spon- 
taneously due to self-grav ity instabilities (see e.g. 
Flohic fc Eracleous ll20Q8l and references therein) 
or could be triggered by close passage of some 
massive object such as another supermassive black 



hole or a star cluster (see e.g. iLewis et al. 1120101 
and references therein). Spiral arms increase flux 
variability of AGN on timescales of a year to sev- 
eral years, but as noted before, they are unable to 
reproduce all aspects of the observed variability. 
However, they are also subjected to fragmentation, 
causing small variations in the flux on timescales 
of several months. The fragments in spiral arms 
can be due to sub-structures in a non-uniform ac- 
cretion disk, such as isolated clumps which could 
pass through the arm and dominate in its emissiv- 
i ty, causing the disc rete "lumps" of excess emission 



( Lewis et al. 1120101 ). The observed variability on 



timescales from few months to several years in the 
difference spectra of some AGN is probably caused 
by such lumps. It is quite possible that some of 
these lumps are long-living and that they do not 
vary significantly in strength, shape, or positio n 
over a period of several years (JLewis et al. II2010I ) . 



As the obtained results show, it is the same case 
with the large bright spots which are responsible 
for amplitude outbursts of the 3C 390.3 H/3 line, 
because they have constant emissivity and they 
are either stationary or spiralling over small dis- 
tances during the period of several years. The only 
feature which significantly varies with time is their 
width. Therefore, these bright spots could be most 
likely explained by the emissivity lumps, caused by 
fragments in spiral arms of the accretion disk. 

5. Conclusions 

We developed a model of the disk perturbing 
region in the form of a single bright spot (or flare) 
by a modification of the power law disk emissiv- 
ity and used this model to simulate the disk line 
profiles. This model has been used to fit the ob- 
served H/3 line of 3C 390.3 observed from 1995 to 
1999. From this investigation we can point out the 
following results: 

1. The model which includes perturbation 
(bright spot) in the accretion disk can suc- 
cessfully explain difference in double peaked 
line profiles, as e.g. higher red peak even 
if we have the standard circular disk. The 
position of a bright spot has a stronger influ- 
ence on one particular part of spectral line 
profiles (such as e.g. its core if the spot is 
in the central part of the disk, or " red" and 
" blue" wings if the spot is located on reced- 
ing and approaching part, respectively). 

2. Using the model for perturbing region we 
were able to successfully model and repro- 
duce the observed variations of the H/3 line 
profile in the case of 3C 390.3, including the 
two large amplitude outbursts observed dur- 
ing the analyzed period. Therefore, the ob- 
served variations of the 3C 390.3 H/3 line 
could be caused by perturbations in the disk 
emissivity. 



3. We found that two outbursts referred 

could 



Shapovalova et al. 



by 
(|200l[ ) could be ex- 
plained by successive occurrences of two dif- 
ferent bright spots on approaching side of 
the disk which are either moving, originat- 
ing in the inner regions of the disk and spi- 
ralling outwards, or stationary. Both bright 
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spots decay by time until they completely 
disappear. 

4. Our results support hypothesis that the per- 
turbations in accretion disk eniissivity are 
probably caused by fragments in the spiral 
arms of the disk. 

The results presented above show that a circu- 
lar disk with perturbations (bright spots) can be 
applied to explain different double peaked line pro- 
files, and can be also used to trace perturbations 
(as well as their characteristics) from the broad 
double peaked line shapes. 
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